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Onset of unsteady horizontal convection in rectangle tank at Pr = 1 
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The horizontal convection within a rectangle tank is numerically simulated. The flow is found to 
be unsteady at high Rayleigh numbers. There is a Hopf bifurcation of Ra from steady solutions to 
periodic solutions, and the critical Rayleigh number Roc is obtained as Roc = 5.5377 x 10® for the 
middle plume forcing at Pr — 1, which is much larger than the formerly obtained value. Besides, 
the unstable perturbations are always generated from the central jet, which implies that the onset 
of instability is due to velocity shear (shear instability) other than thermally dynamics (thermal 
instability). Finally, Paparella and Young's first hypotheses about the destabilization of the flow is 
numerically proved, i.e. the middle plume forcing can lead to a destabilization of the flow. 
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PACS numbers: 47.20.Bp, 44.25. +f, 92.10.af 



Horizontal convection, in which the water is unevenly 
heated at the horizontal surface, was taken as a model 
of abyssal ocean circulation. As the abyssal ocean cir- 
culation plays an important role in climate change, the 
horizontal convection has intensively been explored in re- 
cent years [H, H, Q • It can be set to motion by any small 
temperature gradient, unlike the Rayleigh-Benard con- 
vection. But similar to Rayleigh-Benard convection, the 
horizontal convection may be unsteady at high Rayleigh 
numbers Ra. There is a critical Rayleigh number Rac, 
and the steady flow is unstable and becomes unsteady 
when Ra > Rac. The unsteady flow in horizontal con- 
vection was first found by numerical simulation []| , then 
was observed in the experiment at Ra > 10^^ This 
unsteady flow is proved to be non-turbulent even as 
Ra — > oo, though the flow field seems to be chaotic I*]. 
The investigation on the unsteady horizontal convection 
flow is relatively less, except for [llll, However, they 
have mainly focused on how the turbulent plume main- 
tains a stable stratified circulation. Yet how the horizon- 
tal convection turned to be unsteady remains an elusive 
problem. 

To understand this problem, both Rac for the onset 
of unsteady flow and instability mechanism are of vital. 
Paparella and Young fi\ found Rac « 2 x 10^ at Pr = 1 
in their simulations, which is &|igniflcantly smaller than 
others' results 0, H, 0, B B S Si- example, Rossby 
(1965), Wang and Huang (2005) found the flow is steady 
and stable for Ra < 5 x 10^ in their experiments [1, 
Yet other numerical simulations 0, 0, [1] have not found 
unsteady flows for Ra < 10^. Paparella and Young [l| 
explained the difference of their results from others' as: 

(i) middle plume forcing in the numerical simulations in- 
stead of sidewall plume forcing in the experiments, and 

(ii) lower aspect ratio (H/L = 1/4) in their simulations. 



Both may lead to destabilization of the flow at lower 
Rayleigh numbers. However, their hypotheses have not 
been intensely investigated. According to a recent inves- 
tigation, the flow is still stable for Ra < 10^^ even at 
a much lower aspect ratio {H/L = 1/10) @. Thus, it 
maybe the middle plume forcing that leads to destabi- 
lization at lower Rayleigh numbers. 

Our interest here is to verify their first hypotheses. Is 
the flow with middle plume forcing less stable than the 
sidewall plume forcing? How does the instability occur? 
To investigate these problems, more accurate numerical 
prediction of Rac is need for both forcing cases, for the 
spatial resolution of simulation was very lower in the past 
(e.g. 128 X 32 coarse meshes are used in [l|). Then the 
flow fleld under both middle and sidewall plume forcings 
are compared, which leads to an affirmative answer of 
the above problem. 

Similar to the previous investigations, we consider the 
horizontal convection fiows within the two-dimensional 
domain, and the Boussinesq approximation is assumed 
to be valid for these flows. As shown in FiglU the hor- 
izontal (y) and vertical (z) regimes are < ?/ < L and 
< z < H, respectively. Similar to the depth L 
is taken as reference length scale and A = H/L = 1/4 
denotes the aspect ratio. Taking account of divergence- 
free of velocity field in Boussinesq approximation, the 
Lagrangian streamfunction \1/ and the corresponding vor- 
ticity u! are introduced. The velocity 1? = {v,w), 
where horizontal velocity v = ^ and vertical veloc- 
ity w = respectively. The governing equations in 
vorticity-streamfunction formulation are P, 0, • 
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TABLE I: Comparison of the benchmark solutions from [Tl|] 
and ^T^, Pr — 0.71. 'I'mid, ^'max are the values in the mid- 
point and the maximum of streamfunction, respectively. And 
Nu is average Nusselt number at the heated wall. The reso- 
lution is 80 X 80 meshes for present results. 
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^max 
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\lXi 


10** 


16.386 


16.811 


8.822 


[12] 


10® 


16.386 


16.811 


8.825 


Present 


10® 


16.430 


16.863 


8.828 


mi 


10^ 


29.361 


30.165 


16.523 


[12] 


10^ 


29.356 


30.155 


16.511 


Present 


10^ 


29.605 


30.448 


16.535 



where J(«',(/)) = |f || - H denotes the nonlinear 
advection term. There are two important dimension- 
less parameter in Eq.([T]), i.e. Raylcigh number Ra — 
aT^TgL^ /{kv) and Prandtl number Pr — v/k, where 
g, ax, AT, L, k and v are gravity acceleration, ther- 
mal expansion coefficient, surface temperature difference, 
length of horizontal domain, thermal diffusivity and kine- 
matic viscosity, respectively. Alternatively, Paparella 
and Youngs used vertical length H as length scale, so 
Ra — 64:RaH, where Ran is the vertical Rayleigh num- 
ber by using vertical length H as unit 

More specifically, we consider the horizontal convection 
in a rectangle tank at Pr = 1. The tank has same veloc- 
ity boundary condition as that in i.e. free slip and no 
shear stress at the walls. In addition, two different sur- 
face forcings are used, which are central symmetric. One 
is middle plume forcing as T = [1 -|- cos(27ry)]/2 the 
other is sidewall plume forcing as T = [1 — cos(27r?/)]/2 
Comparing these with one cell forcing T = 
cos(7ry/2) there are two symmetric cells in the flow 
field under such forcings (e.g. Fig[T]and Fig[6]), when the 
flow is symmetrically steady and stable. In addition, the 
middle plume forcing in the left cell is the same with the 
sidewall plume forcing in the right cell (see FigUbehind). 
Thus in the steady flows, both forcings will lead to the 
same flow patterns except for a position shift, which is 
proved by the following investigation. 

There are two important quantity describing the cir- 
culation, i.e. the non-dimensional streamfunction max- 
imum and the non-dimensional heat flux. The non- 
dimensional streamfunction maximum ^'max = ^max/'^i 
where ^'J^ax the maximum of the dimensional stream- 
function. 

The above Eq. ([T|) is solved with flnite different method 
in non-uniform grids. Crank-Nicholson scheme and 
Arakawa scheme [e.g. [l^ [ij] are applied to discretize 
the linear and nonlinear terms, respectively. Comparing 
to the other schemes, Arakawa scheme is more accurate 
but more expensive, and it has also been applied to hor- 
izontal convection flows at high Rayleigh number ^ Q . 
Table [T] shows the validation of the scheme with nature 
convection problem. A flne spatial resolution mesh of 
512 X 128 is used to eliminate numerical instability. 
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FIG. 1: The flow stream function (a) and temperature field 
(b) of Ra — 5 X 10*. It is steady and stable and symmet- 
ric with middle plume forcing, solid and dashed curves for 
positive and negative values, respectively. 



First, the middle plume forcing is considered, which is 
steady and stable for Ra < 5.5 x 10*. Fig[l] shows the 
flow field (a) and temperature field (b) of Ra = 5 x 10* 
with ^'max = 59.83, in which the flow is symmetric, 
steady and stable. In this case, the center line symmet- 
rically separates the flow flcld into two parts, like a free 
slip wall. There is a vigorous downward jet in the cen- 
ter of tank corresponding to the middle plume forcing 
(Figd^), where the vertical velocity field has a minimum 
oiw = —2513 (Figl^b). The center jet leads to the clock- 
wise and anticlockwise plume cells in the left and right 
part of tank, respectively. In the left circulation cell, 
the flow sinks quickly along the center line and upwells 
clockwise along the left side wall with relatively slower 
speed, which can be also seen from the vertical velocity 
of the flow (Fig[2]). Besides, there are two horizontal jets 
respectively near top and bottom walls in the left circula- 
tion cell (Figl2^). Totally, there are 2 horizontal jets near 
wall and a vertical jet at the center in each cell. Contract 
to the flow fleld, the temperature field is very simple. An 
obvious boundary layer exists near the surface in temper- 
ature field, which leads to a 1/5-power law of Ra for heat 
flux [e.g. IB, S And below the temperature bound- 
ary layer, the temperature is almost homogeneous due to 
the convection. Thus there is a very strong stratification 
near the surface {dT/dz ~ Ra'^/^) but a very weak strat- 
ification in other region {dT/dz ~ 0). As the above case 
is stable, so the critical Rayleigh number must be larger 
than 5 x 10*, which is significantly larger than the value 
obtained before 

To find the critical Rayleigh number i?ac, the growth 
rate of perturbation (pit) is calculated numerically. And 
(f){t) is assumed to satisfy (t>{t) — e°'*(/)(0), where cr = (7^ + 
iui is the complex growth rate of disturbance. It is found 
that the onset of unsteady flow is at Rac = 5.5377 x 10*, 
as shown in FigO For Ra — 5.53 x 10*, the flow is stable 
and the growth rate is approximately Ur = —0.12. But 
the flow is unstable and the growth rate is approximately 
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FIG. 2: The horizontal (a) and vertical (b) velocity fields of 
Ra = 5 X 10*. It is steady and stable and symmetric with 
middle plume forcing, solid and dashed curves for positive 
and negative values, respectively. 




FIG. 3: Growth rate ar (solid) and ai (dashed) vs. Ra, re- 
spectively. 



ar = 0.03 for Ra = 5.54 x 10*^. Thus the critical Rayleigh 
number Ruc is obtained 5.53 x 10^ < Rqc < 5.54 x 10*. 
The accurate value of Rqc — 5.377 x 10® is obtained by 
interpolating from the above result. Moreover, the onset 
of unsteady flow is found to occur via Hopf bifurcation. 
As Figl3] shows, the image part of growth rate is nonzero 
and the eigenmode of perturbation is periodic. This Hopf 
bifurcation of the horizontal convection has not been re- 
ported yet, and previous investigations dealt only with 
chaotic flows. 

Meanwhile, the evolution of the perturbation vorticity 
fields during the first half period at t — {sl), t = T/8 
{h),t = T/4: (c) and t = 3T/8 of Ra = 5.54 x 10® are de- 
picted in Figdl respectively. The perturbation vorticity 
fields are symmetric about centerline, which implies that 
the horizontal velocity is nonzero at centerline. It can 
be seen that the perturbation tripole A (the shadowed 
ellipse in FigH^) is generated from central downward jet, 
then propagates and amplifies along the central jet down- 
ward to the bottom wall (Fig|3}D,c,d). When tripole A ap- 
proaching to the bottom, it becomes weaker and weaker 
and breaks into two parts: the left and the right near the 




FIG. 4: The perturbational vorticity fields at t = (a), 
t = T/8 (b), t = T/4 (c) and t = 3T/8 of Ra = 5.54 x 10^ 
solid and dashed curves for positive and negative values, re- 
spectively. 
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FIG. 5: The vorticity of Ra — 5.54 x 10* with vertical veloc- 
ity w (shadowed as u; > 1800), solid and dashed curves for 
positive and negative values, respectively. 



bottom, which can be seen from the evolution of tripole 
B (the shadowed rectangle in FigllK). And the mean flow 
advects the broken vortexs horizontally along the bottom 
wall (Fig[3]D,c,d). Then in the second half period, a re- 
verse tripole will generate the same place of vortexes A 
at t = T/2, and the same story repeats for it, which is 
omitted here. In short, the perturbations generate and 
amplify in the central vertical jet, but are propagated 
and weakened along the horizontal wall. 

Further investigation shows that the instability of flow 
occurs due to shear. First, as we noting, the instabilities 
always occur in the center and propagate along the mean 
flow. Second, this trigger place locates in the area where 
there is a vigorous jet with strong shear (see e.g. Figl^b 
and Fig[5]) (FiglS]). As the stratification is very weak here 
(see e.g. FigUlD), the flow in this region is dominated by 
momentum dynamics other than thermal dynamics. All 
these imply that the onset of instability leading to un- 
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FIG. 6: The flow field (a) and temperature field (b) of Ra — 
5 X 10*, which are steady and stable and symmetric with 
sidewall plume forcing. 



steady flow is due to shear instability at larger Rayleigh 
numbers, which is much different from Rayleigh-Benard 
instability. However, shear is not the sufficient condition 
for instability. For example, the instability near the top 
surface is suppressed due to strong stratification (FiglTJa), 
though both the velocity (Fig[2^) and the shear (FiglS]) 
near top surface are still very large. In words, the onset 
of instability is due to velocity shear (shear instability) 
other than thermally dynamics (thermal instability). 

Now, the sidewall plume forcing is considered. Fig [5] 
shows the flow field and temperature field of Ra = 5x10^, 
in which the flow is symmetric, steady and stable like 
that under middle plume forcing. There are two vigor- 
ous downward jets near the walls corresponding to the 
sidewall plume forcing (FiglH^). As mentioned above, 
the sidewall plume forcing will lead to exactly the same 
flow pattern as the middle plume forcing does except for 
a position shift, which can be seen from FigUand FiglS] 
As the flow is stable, the center line like a free slip wall 
symmetrically separates the two cells. The left cell in 
FiglT] is exactly the same with the right cell in FigEl 
However, the flow is much more stable with the sidewall 



plume forcing. And the critical Rayleigh number is found 
Ruc 1.85x10^° (with 768x192 meshes) in this case. As 
mentioned above, the flow is much more stable with the 
sidewall plume forcing than that with the middle plume 
forcing, though both forcings lead to the same flow pat- 
terns. This is very interesting, and can be understood 
from the mechanism of instability. 

It's found that the rigid wall suppresses the perturba- 
tion, which leads a more stable flow with the sidewall 
plume forcing than that with the middle plume forcing. 
As the loss of stability is due to strong velocity shear in 
the center in horizontal convection, the smaller the shear 
is, the more stable the flow is. In the case of middle 
plume forcing, the perturbation with nonzero horizontal 
velocity occurs at the most vigorous downward jet. And 
the perturbed flows cross the center line and propagate 
downstream. However, in the case of sidewall plume forc- 
ing, these crossing flows are suppressed by rigid walls. 
So that the critical Rayleigh number is much larger in 
this case. Paparella and Young (2002) hypothesized that 
middle plume forcing may lead to a destabilization of the 
flow. Here this hypotheses is proved both physically and 
numerically. 

In conclusion, the onset of unsteady flow is found to oc- 
cur via a Hopf bifurcation in the regime of Ra > Ra^ = 
5.5377 X 10® for the middle plume forcing at Pr — 1, 
which is much larger than the previously obtained value. 
Besides, the onset of unsteady flow is due to shear in- 
stability of central downward jet. Finally, the flrst hy- 
potheses of Paparella and Young (2002) for instability is 
numerically approved, i.e. the middle plume forcing can 
lead to a destabilization of the flow at relatively lower 
Rayleigh numbers. 
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